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Recently it has been shown [I] that fluorescence telescopes with a large field of view can indirectly 
probe the scale of supersymmetry breaking. Here we show that depending on their ability to fight 
a large background, multi-Km 3 volume neutrino telescopes might independently probe a similar 
breaking scale region, which lies between ~ 10 s and ~ 5 x 10 6 GeV. The scenarios we consider 
have the gravitino as the lightest supersymmetric particle, and the next to lightest (NLSP) is a long 
lived slepton. Indirect probes complement a proposal [2] that demonstrates that 1 Km telescopes 
can directly probe this breaking scale. A high energy flux of neutrinos might interact in the Earth 
producing NLSPs which decay into taus. We estimate the rate of taus, taking into account the 
regeneration process, and the rate of secondary muons, which are produced in tau decays, in multi- 
km 3 detectors. 

PACS numbers: 14.80.Ly, 12.60. Jv,95.30.Cq 



I. INTRODUCTION 

Extensions of the standard model of particle physics 
(SM) are currently being probed at the Large Hadron 
Collider (LHC). At the same time, features of these ex- 
tensions can be probed by non accelerator experiments. 
It has recently been shown [1] that the supersymmetry 
breaking scale (yF) can be probed indirectly by large 
fluorescence telescopes, such as the JEM-EUSO Obser- 
vatory. This scale, as well as Universal Extra Dimensions 
scenarios [3], can also be directly probed by 1 Km 3 vol- 
ume neutrino telescopes 0|4], such as IceCube. Here we 
show that depending on their ability to fight a large back- 
ground, multi-Km 3 volume neutrino telescopes might 
also indirectly probe the scale of supersymmetry break- 
ing. 

We investigate scenarios where supersymmetry (susy) 
is broken at \f~F~ < 10 10 GeV , and the lightest supersym- 
metric particle is the gravitino. In many of these scenar- 
ios, including Gauge Mediation Susy Breaking models [5], 
the NLSP is a right-handed stau, whose lifetime is [3J: 

where nif R is the stau mass. The NLSP range depends 
on y/F, and for yF > 5 x 10 6 , they travel a long way 
before decaying while for lower values they decay after a 
short travel. It has been shown [2, 4 J that neutrino tele- 
scopes can directly probe this scale, probing scenarios 
where 5 x 10 6 < y/F < 5 x 10 s GeV, while large fluo- 
rescence telescopes [DEj, such as the JEM-EUSO Obser- 
vatory [7], can probe scenarios where NLSPs decay and 
VF < 5 x 10 6 GeV. 
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In this article we show that multi-Km 3 volume neu- 
trino telescopes can probe the same vF region as large 
fluorescence telescopes. NLSP decays will yield a few 
events per year in the projected KM3Net neutrino tele- 
scope [8] , or in an expanded version of IceCube [9] . How- 
ever, while fluorescence telescopes have a clean signature 
of these events, probes by neutrino telescopes will de- 
pend on the ability to discriminate taus (rs) produced in 
NLSP decays from a large background. 

A detectable rate of NLSPs in a 1 Km 3 neutrino tele- 
scope, such a IceCube, can be produced by high energy 
neutrinos interactions in the Earth. Here we consider the 
same NLSP production mechanism, but analyze scenar- 
ios where these decay, probing a complementary breaking 
scale, where \[F < 5x 10 6 GeV. 

High energy ts are produced in NLSP decays and un- 
dergo a regeneration process, where r neutrinos (y T ) pro- 
duced from r decay, charge current interact in the Earth 
and produce a new r. This process will expand the frac- 
tion of the Earth volume within which the r can decay 
and still reach the detector. We also consider the flux of 
muons produced when these secondary ts decay. 

The caveat in probing y/F [5J comes from the fact that 
NLSPs lose less energy and have a much larger range 
while transversing the Earth when compared to SM lep- 
tons. This larger range compensates its lower production 
cross section. Here, when NLSP decays are considered, 
an extra decaying volume is gained due to the r regener- 
ation process. 

Multi-km 3 neutrino telescopes proposals exist, either 
as a new telescope in the Mediterranean Sea [5] or as an 
expansion of IceCube [9] . Here we show that depending 
on their final design and their ability to fight a large 
background, they might be able to probe the y/F region 
which allows NLSPs decays in the Earth. 

In the next section we describe our simulation of NLSP 
production, propagation and energy loss, which repro- 
duces the analysis described in [2j|4]. Follows the descrip- 
tion of our simulation of NLSP decays with r production; 
r decay and regeneration process as well as muon pro- 
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duction, which is the same used in PQ. The event rate 
in neutrino telescopes will depend on its design. In Sec- 
tion |III| we discuss this point, and determine the rates 
of t hese events in a 1 Km 3 neutrino telescope. In sec- 
tion IV the rate is determined for multi-km 3 telescopes. 
In Section |"IV A| we analyze r specific signatures, as well 
as a cosmological background. Finally we present our 
conclusions. 



II. NLSP PRODUCTION 

As shown in [2], NLSPs can be produced from the in- 
teraction of high energy neutrinos in the Earth. Although 
the production cross section related to this process is even 
lower than the one for SM leptons, the flux of NLSPs will 
still be considerable due to its much larger range, as de- 
scribed in the previous section. 

Here we take the Waxman-Bahcall (WB) [TO] upper 
limit on the neutrino flux hitting the Earth as a reference 
for our analysis. Our result is normalized by this upper 
limit, and can be translated to any other flux if properly 
rescaled. We consider neutrino production from optically 
thin sources, and the WB limit (for a maximized cosmo- 
logical evolution of sources) is given by 



d<t> v 
dE„ 
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where E v is the neutrino energy. It is important to note 
that although this limit was set for muon and electron 
neutrinos, our analysis is independent of the initial neu- 
trino flavor. Any flavor of the left-handed sleptons, pro- 
duced from neutrino interactions, immediately decay pro- 
ducing a chain that will always end up with right-handed 
staus (the NLSPs). For this reason, neither the initial 
neutrino flavor nor cosmogenic mixing due to oscillations 
alter our results. 

NLSPs will be produced by the interaction of cosmo- 
genic neutrinos in the Earth and the production cross 
section is determined as in 4 . It follows the same pat- 
tern as for SM lepton production by charge current inter- 
actions, where now the up or down type quark is a squark 
(q) and the chargino is the mediator in the t-channel. We 
also account for the sub-dominant neutralino exchange. 
From this interaction a left-handed slepton (Zj,) and a q 
are produced and promptly decay, always producing two 
NLSPs at the end of the decay chain. The production en- 
ergy threshold depends on rriq and mj- and is ~ 10 5 GeV, 
and the NLSP initial energy is ~ E u /6 g] . These NLSPs 
are typically right handed staus (tr) which travel in par- 
allel through the Earth. We take the chargino mass as 
ma = 250 GeV, and m f = 150 GeV, m fR = 250 GeV for 
the left-handed stau and the NLSP respectively, and con- 
sider three possibilities for the squark mass, rriq = 300, 
600 or 900 GeV. There are constrains to the tr from 
big-bang nucleosynthesis [11] . 

NLSP production depends on the probability that neu- 
trinos will interact in the Earth, and we use the Earth 
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FIG. 1. Number of NLSP events per year in a km 3 neutrino 
telescope versus the scale of susy breaking yF. It is assumed 
an initial neutrino flux equivalent to the WB limit. The curves 



are shown for different rriq values. 



density profile model described in [TU H3] • Our simula- 
tion of NLSP propagation includes its energy loss, due 
both to ionization and radiative processes. It follows the 
analysis described in [U [2] , where it is shown that the 
NLSP energy degradation due to radiative losses is sup- 
pressed when compared to SM leptons. 

We check our NLSP production and propagation sim- 
ulation by reproducing the direct detection rate and the 
energy distribution shown in [21 H] . 



NLSP DECAY AND EVENT RATES IN 1 
KM 3 NEUTRINO TELESCOPES 



III. 



Here we consider scenarios where v F is such that NL- 
SPs decay after a short travel through the Earth. This 
distance equals 7cr, and can be determined from Equa- 
tion [l] They decay into a r and a gravitino, where the r 
is much heavier than the gravitino. Figure [T] shows the 
number of NLSPs per year in a km 3 neutrino telescope 
(positioned ~ 2 Km deep in the Earth) versus y/F, and 
it is clear that NLSPs decay for y/F < 5 x 10 6 GeV. 

In our Monte Carlo simulation ?9(10 5 ) neutrinos 
isotropically distributed reach the Earth. These are gen- 
erated with energies between 10 5 to 10 12 GeV, and are 
normalized by the WB limit. NLSPs are then produced 
according to the neutrino interaction probability distri- 
bution, which depends on the Earth density profile and 
the NLSP production cross section, and propagated, tak- 
ing into account their energy loss, as described in the 
previous section. The neutrino interaction probability is 
convoluted by its survival probability to account for SM 
lepton production. 
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FIG. 2. Average number of r regenerations. Only the rs 
which reach a 1 km 3 detector are taken into account. Results 
are shown for three values of rriq, and as expected show no 
dependency on this parameter. 



Once the NLSP production and propagation is sim- 
ulated, NLSPs are randomly selected to decay accord- 
ing to their decay probability distribution, Pd = 1 — 
exp(—mx/Ef L cT). The decay generate rs isotropically 
distributed in the center of mass frame, which are 
boosted to the laboratory frame, mainly following the 
same direction as the NLSPs. rs will then decay and 
generate v T s and a regeneration process might take place. 
In 18% of the r decay, a muon will be produced and, if 
close enough to the detector, its propagation will also be 
simulated. 

The regeneration process depends on the r energy. The 
fraction of energy carried by the v T is determined as de- 
scribed in [151 116) . The average number of regenerations 
are shown in Figure [2] as a function of the neutrino en- 
ergy E V1 for rs that reach a 1 Km 3 detector positioned 
as IceCube. 

We note that two NLSPs are always produced from 
the neutrino interaction in the Earth. These NLSPs will 
decay at different points in the Earth and the generated 
rs will not follow the same coincidence pattern in the 
detector as the NLSPs. 

In order to understand the possibility of NLSP indi- 
rect detection by neutrino telescopes, the detector shape 
is relevant. The r average path and number of regen- 
erations depend strongly on its energy. Up to 10 6 GeV 
the average path length is much smaller than a Km, but 
larger energies imply that this distance might be longer 
than the detector size. Events arriving horizontally with 
respect to the detector can go through fewer regeneration 
processes than the ones coming from below. Therefore, 
it is important to determine the event rate as a function 
of the detector volume as well as to take into account the 



FIG. 3. Energy distribution of produced f pairs (blue square), 
the fraction of these that reach the detector (red dots) assum- 
ing y/F = 10 7 GeV and of rs (pink triangles) and fj,s (green 
triangles) as well as their sum (black stars) that reach the de- 
tector assuming \/~F = 10 5 GeV. Here we take rriq = 300 GeV. 
The distributions for rriq = 600 and 900 GeV, are very similar 
in shape with lower number of events. 



detector shape. One important point is that due to the 
regeneration process, the effective volume for NLSP pro- 
duction will be enhanced when compared to direct NLSP 
detection. 

We consider a 1 km 3 cylindrical detector, with a 
0.564 Km radius, 1 Km height, and buried 2 Km deep in 
the Earth. Our results are therefore relevant both for Ice- 
Cube [S] and Km3Net [5] neutrino telescopes. We show 
our results as a function of two arbitrary values of vF, 
10 5 and 10 7 GeV, where the lower value represents the 
maximum NLSP decay rate and the upper one the region 
where NLSPs are stable while transversing the Earth. 

Figure [3] shows the energy distribution for events reach- 
ing a 1 Km 3 detector. Up to neutrino energies of 10 7 GeV 
the r rate is smaller than the fi one and this relation 
inverts at ~ 10 8 GeV. At lower energies the r decays 
almost instantly, which lower it's detection probability 
when compared to larger energies where it travels for a 
while before decaying. For E v > 5 x 10 9 GeV all rs will 
reach the detector. This figure assumes rriq = 300 GeV, 
whereas the distributions for rriq = 600 and 900 GeV, are 
very similar in shape with lower number of events. The 
total rate of events in a Km 3 detector per year, for the 
three values of rriq, are shown in Table ^ It is important 
to note that the number of events in larger detectors will 
not necessarily scale linearly with these results, as we will 
show below, due to the importance of the shape of the 
detector. 

The r and [i energy at the detector is shown in Fig- 
ure |U It can be seen that most rs that reach the detector 
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mg (GeV) 


V (km 3 ) 


r (y- 1 ) 


t (y- 1 ) 


a* (y 1 ) 


V + t (y -1 ) 


300 


1 


27.7 


0.07 


0.21 


0.28 


300 


21 


521.5 


1.60 


2.36 


3.95 


300 


65 


1450.5 


4.41 


5.85 


10.26 


600 
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4.5 


0.02 


0.04 


0.06 


600 


21 


100.6 


0.51 


0.47 


0.98 


600 


65 


260.0 


1.47 


1.27 


2.74 


900 


1 


1.1 


0.01 


0.01 


0.02 


900 


21 


28.2 


0.22 


0.15 


0.37 


900 


65 


73.0 


0.65 


0.41 


1.05 



TABLE I. t pair, r and u, rates per year in a Km neutrino 
telescope at ~ 2 Km deep, for neutrino telescopes of 1, 21 
and 65 Km 3 volumes. Results are given for 3 rriq, assuming 
y/F = 10 r (10 5 ) GeV for f pairs (rs and fis). 



V = 1 km 3 , (F = 1x10 5 GeV 




9 10 11 
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FIG. 4. Distribution of r and /i energy at a 1 Km 3 detector for 
VF = 10 5 GeV. Results are shown for m q = 300 GeV. Other 
values of rriqS have similar distributions with less events. 



have larger energies than the /its, reflecting the distribu- 
tion in Figure [3] 
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FIG. 5. Same as Figure [3] but now for extended neutrino 
telescopes. Top plot shows results for 21 Km 3 and bottom 
for 65 Km 3 volume telescopes. Both show results for rtiq = 
300 GeV. It also shows the cosmological muon background 
(black * line) energy distribution. 



IV. DECAY RATES IN MULTI-KM 3 
TELESCOPES 

As can be seen in Table |TJ the indirect probe of the susy 
breaking scale is not promising at 1 Km 3 neutrino tele- 
scopes. However there are studies [§] and proposals [5] 
of multi-Km 3 detectors. Here we determine the r and fi 
rates from NLSP decays in these detectors and analyze 
how well they can probe this scale. 

We take a 2 and 4 Km radius extension to our current 
cylindrical configuration (as proposed in [9j for an Ice- 
Cube extension) as our case scenarios. This leads to 21 
and 65 Km 3 volumes. It is important to note that due 



to the regeneration process, the increase in the indirect 
detection rate will depend on the shape of the extended 
detector. 

In Figure [5] we show the r and /i rates in such tele- 
scopes, where we assume y/~F — 10 5 GeV, representing 
scenarios where all NLSPs decay. The event rate per year 
for these extended volumes are shown in Table HI The dis- 
tribution for the r and /x arrival energy at the detector 
for these events are very similar to the ones shown in 
Figure [4] 

As shown in Table [I] the r and /i rates in extended 
neutrino telescopes are significant, with a few events per 
year for rriq — 300 GeV and 1 event per year in the 
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larger volume for Trig = 900 GeV. These rates are at the 
same level as NLSP direct detection [2] [4] and indirect 
dark matter detection [T7] in 1 Km 3 neutrino telescopes. 
Therefore one needs to determine possible backgrounds 
for these events. 



A. Backgrounds in Multi-Km 3 Telescopes 

Under the scenarios we are considering, while NLSP 
direct detection has a clean and very distinctive signa- 
ture, with two NLSPs transversing the detector simul- 
taneously, NLSP indirect detection has a considerable 
background. It comes both from atmospheric as well as 
from cosmological neutrinos, where the latter is also the 
source of NLSPs. In the indirect case, each produced 
NLSP will decay at different times and therefore almost 
all rs and /is will transverse the detector as a single event. 

The number of events in detectors of different volumes 
is shown in Table [IJ Figure [5] shows the energy distri- 
bution of cosmological /is per year in both 21 and 65 
Km 3 telescopes. Muons from atmospheric neutrinos have 
lower energy compared to the cosmological neutrinos and 
can be significantly reduced. However cosmological neu- 
trinos will yield an almost unreducible background. The 
only way to distinguish /is coming from NLSP decays, 
would be through a time correlation with the NLSP pro- 
duced rs. 

However there is also a small but significant number 
of rs for extended telescopes. In the best case scenario, 
four taus per year can be seen in the 65 Km 3 volume 
telescope. There are many techniques to distinguish rs 
from /is in neutrino telescopes, as described in the next 
section. Most of these signatures are favored by an in- 
creased telescope size. As an example, in the case of the 
double-bang [18] signature, two cascades can occur inside 
the detector. For a 1 km 3 telescope there is an tight r en- 
ergy window, in order to contain these cascades, while a 
larger detector has a much larger efficiency for this signa- 
ture. There is an ongoing effort to improve discrimination 
between /is and rs in neutrino telescopes |19|. Once this 
discrimination is done efficiently, the background for rs 
generated by NLSP decays will depend on neutrino oscil- 
lations. If one assumes that cosmological neutrinos oscil- 
lates, always arriving in a 1 : 1 : 1 ratio of v ei v^ and v T s, 
there will also exist an almost irreducible r background 
(see next section). In this case, /is and rs produced by 
NLSP decays have to be time correlated in order to re- 
duce the background. Otherwise, NLSP produced rs will 
be background free. 



V. t SIGNATURES IN NEUTRINO 
TELESCOPES 

In order to check our regeneration simulation, we gen- 
erated rs from a distribution of 10 6 ^ r s isotropically hit- 
ting the Earth. The regeneration process was then sim- 
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FIG. 6. Fraction of generated events which produce a specific 
r signature as labeled. In red is shown the total number of 
rs that reach the detector. A 1 Km 3 volume is assumed. 



ulated as described in Section III We analyzed the fol- 
lowing r signatures [IS] [2D] : double bang, where the v T 
charge current interacts in the detector producing a r 
plus a hadronic shower and the produced r subsequently 
decays producing a second shower; lollipop, where a r 
enters the detector and decays inside it producing one 
shower; inverted lollipop, where one shower is produced 
together with a r which in its turn goes through the de- 
tector, and finally the tautsi pop signature, where as in 
the double bang, two showers are produced but too close 
to each other looking like a single shower. 

Figure [6] shows the fraction of generated events corre- 
sponding to each r signature assuming a 1 km 3 volume, 
and the number of rs that go straight through the detec- 
tor. Table [IT] shows the fraction of events for each specific 
signature, integrated in energy and the number of events 
assuming that 1 /3 of the WB upper limit corresponds to 
v T s. Both cases integrates the energy above 2 x 10 6 GeV, 
corresponding to the minimal energy with which most rs 
arrive at the detector (see Figure [4| . These results are 
consistent with [2"T] . 

Table [fT] shows that if 1/3 of the cosmological neutrino 
flux arrives at the Earth as there will also be a signif- 
icant cosmological background for the NLSP generated 
rs. The only chance for extended neutrino telescopes to 
probe decaying NLSPs will therefore be a time correla- 
tion between both rs or a r and a /i produced from NLSP 
decay. 



VI. CONCLUSIONS 

We have shown that multi-Km 3 neutrino telescopes 
are potentially sensitive to indirectly probe the scale of 
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TABLE II. Fraction of generated events with specific r sig- 
nature (see text), for 1, 21 and 65 Km 3 neutrino telescopes; 
and integrated number of events assuming an initial v r flux 
equivalent to 1/3 of the WB limit. Only events with energies 
above 2 x 10 6 were accepted, since most NLSP produced rs 
arrive the detector with energies above this value. 



V = 1 km 3 


Limit 


Total 


Lollipop 


Inv 
Lollipop 


Tautsipop 


Straight 
Through 


Double 
Bang 


/«(io- IU ) 


7.9 


0.37 


0.56 


0.009 


4.9 


0.3 


WB 


31 


1.5 


2.2 


0.037 


19.5 


1.3 


V = 21 km 3 


/«(io- lu ) 


270 


29 


36 


0.3 


120 


40 


WB 


1052 


114 


143 


1.3 


457 


158 


V = 65 km 3 


/«(10- 10 ) 


860 


83 


130 


1.2 


300 


200 


WB 


3420 


328 


517 


4.7 


1180 


790 



susy breaking. However although a significant number 



of NLSP generated rs will trigger these telescopes, this 
probe depends on the ability to distinguish these rs from 
a large background. This background is composed of rs 
produced by cosmological neutrinos, that will be signifi- 
cant if one assumes an equal neutrino flavor ratio at the 
Earth, due to oscillations. If for any reason oscillations do 
not produce tau neutrinos, extended neutrino telescopes 
are able to probe the susy breaking scale in the range 
10 5 H vF < 10 6 GeV. Large fluorescence telescopes can 
also indirectly probe [1] the same y~F range but with the 
advantage of a background free signature. Both these 
proposed indirect measurements complement direct susy 
breaking probes that can be done by 1 Km 3 neutrino 
telescopes [HE]- 
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